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ABSTRACT 

Halo-based models have been successful in predicting the clustering of matter. How¬ 
ever, the validity of the postulate that the clustering is fully determined by matter 
inside haloes remains largely untested, and it is not clear a priori whether non-virialised 
matter might contribute significantly to the non-linear clustering signal. Here, we in¬ 
vestigate the contribution of haloes to the matter power spectrum as a function of 
both scale and halo mass by combining a set of cosmological N-body simulations 
to calculate the contributions of different spherical overdensity regions, Friends-of- 
Friends (FoF) groups and matter outside haloes to the power spectrum. We find that 
matter inside spherical overdensity regions of size i? 200 ,mean cannot account for all 
power for 1 < k < lOOftMpc - , regardless of the minimum halo mass. At most, it 
accounts for 95% of the power (k > 20/iMpc -1 ). For 2 < k < 10/iMpcU 1 , haloes with 
mass M 2 oo,mean /S 10 11 ft -1 Mq contribute negligibly to the power spectrum, and our 
results appear to be converged with decreasing halo mass. When haloes are taken 
to be regions of size f?200,crit, the amount of power unaccounted for is larger on all 
scales. Accounting also for matter inside FoF groups but outside -R 2 oo,mean increases 
the contribution of halo matter on most scales probed here by 5 — 15%. Matter in¬ 
side FoF groups with M 2 oo.mean > 10 9 ft -1 Mq accounts for essentially all power for 
3 < k < 100/iMpcU 1 . We therefore expect halo models that ignore the contribution 
of matter outside i? 2 oo,mean to overestimate the contribution of haloes of any mass to 
the power on small scales (k > 1 ft, Mpc -1 ). 

Key words: cosmology: theory, large-scale structure of Universe - galaxies: haloes 


1 INTRODUCTION 

The matter power spectrum, a measure of how matter 
clusters as a function of scale, is a key observable of 
our Universe. As future weak lensing experiments which 
aim to measure this quantity with unprecedented accu¬ 
racy, such as DElfl LSS'10. Euclic0 and WFIRST0, draw 
closer, the precision with which the theoretical matter 
power spectrum is predicted must also increase. Currently, 
some of the largest uncertainties on fully non-linear scales 
come from our incomplete understanding of galaxy forma- 
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tion (e.g. [van Daalen et al. !H), which can cause large 
unwanted biases in the cosmological parameters derived 
from observations. We may be able to correct for these 
biases using independent measurements of, for example, 
the large-scale gas distribution, and/or by marginalising 
over these uncertainties using a halo model based ap¬ 
proach. However, for the largest of the future surveys more 


than currently exist will be needed le.E. Semboloni et al. 

2011 


Semboloni. Hoekstra & Schave 2013!, Zentner et al. 

2013 

Nataraian et al. 20141, Mohammed & Seliak 20141. 


Eifler et al.l 20141 '). 


But even assuming that we can somehow account 
for the effects of galaxy formation on the distribution 
of matter, significant challenges remain before we are 
able to predict the matter power spectrum with the 
sub-percent accuracy needed to fully exploit future mea- 
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surem ents llHuterer fe Ta kada 200l|, iHearin. Zentner fe Mai 
l2012h . These include converging on the “true” simulation 
parameters in N-bod y codes, althou gh these too can be 
marginalised over JSmith et all 120141 ). However, with each 
such marginalisation one should expect the constraining 
power of the observations to be reduced. 

Direct simulations are not the only way to obtain the¬ 
oretical predictions for the matter power spec trum . Other 
avenues, such as the an alyti cal halo model (e. g. Seliakll2000l . 
IPeacock fc Smithl l2000l : see ICoorav fc Shet'hl 120021 for a re¬ 
view), exist and are widely used in clustering studies. The 
halo model is based on the assumption that all matter is 
partitioned over dark m atter haloes, which finds its origin 
in the model proposed bv Press_^Schechterl d 19T4I . hereafter 


PS), later extended bv lBond et al. 1 199lh . The PS formal¬ 


ism is based on the ansatz that the fraction of mass in haloes 
of mass M(R) is related to the fraction of the volume that 
contains matter fluctuations 5r > 5 C rit, where <5 is the linear 
density contrast, R is the smoothing scale and <5 C rit is the 
critical, linear density contrast for spherical collapse. If the 
initial field of matter fluctuations is known, then a halo mass 
function can be derived from this ansatz, which together 
with a model for the bias b(M) (the clustering strength of 
a halo of mass M relative to the clustering of matter) and 
a description of halo density profiles fully determines the 
clustering of matter. 

Much work has been done to improve the predic¬ 
tions of the halo model since its introduction. More 
accurate mass functions have been derived based on, 


200111, fits to N-body simulations (e.g. 

Jenkins et al. 200ll. 

Tinker et al. 2008. 

Bhattacharva et al. 

201 ll. Angulo et al. 

2012 

, Watson et al. 

20131; see Murray, Power & Robotham 

2012 

for a comparison of different 

models) and sim- 


ulatio ns taking into accoun t the effects of baryons 


(e.g. IStanek. Rudd fe Ev rardl 12 0091, ISaw ala et all 


Martizzi et al. 


Velliscig et al. 


I20l4 Cui. Borgani fe Murantel [201 4 
120141 ). Similarly, much effort has gone 


2013 


tions (e.g. Sheth & Tormen 

1999. Seljak & Warren 

2004 

Smith. Scoccimarro & Sheth 

[2007], 

Reed et al. 


2009 

Pillepich, Porciani & Hahn 20ld, Tinker et al. 

201C 


and concentration-mass relations for halo profiles (e. 


Bullocl^eMd. 2001 Neto et all I20071. Duffy et al.l 200 


JYI. IJJulty . 

Maccio . Dutton & van den Bosch 20081 ILudlow et al.l 
20141 '). Current halo models may incorporate addi¬ 
tional ingredients like triaxiality, substructure, halo 
exclusion, p rimordial non -G auss ia nity and baryonic ef¬ 
fects 


Smith fe Watts 


_ (e.g. Isheth fe Jainl 200 31. _ __ 

Giocoli et al.l hold. ISmith. Desiacaues fe Marian 


] 

2005 


2011 


Gil-Marfn. Jimenez & Verde 1201 ill . and fitting formulae 
based on the ha l o model have also be en developed (e.g. 
ISmith et al.ll2003l . iTakahashi et al.ll2012l ) . 

However, the validity of the postulate that the cluster¬ 
ing of matter is fully determined by matter in haloes re¬ 
mains relatively untested. Even though matter is known to 
occupy non-virialised regions such as filaments, their mass 
may simply be made up of very small haloes itself, although 
recent results indicate that part of the dark matter ac¬ 


creted onto_Jialoes_Js genuinely_smooth ( Angulo_&_Whitc 


201C , iFakhouri fe Mai l20ld . iGenel et all 120101 . Wang et al 


20111 ). Either way, it is not clear a priori whether this non- 


Table 1 . The different simulations employed in this paper. From 
left to right, the columns list their name, box size, particle mass 
and maximum proper softening length. All simulations were run 
with only dark matter particles and a WMAP7 cosmology. 


Name 

Box size 
[, h -1 Mpc] 

Particle 

number 

m dm 

[h- 1 M 0 ] 

%ax 
[/i _1 kpc] 

L 400 

400 

1024 3 

4.50 x 10 9 

4.0 

L200 

200 

1024 3 

5.62 x 10 8 

2.0 

L100 

100 

512 3 

5.62 x 10 8 

2.0 

L050 

50 

512 3 

7.03 x 10 7 

1.0 

L025 

25 

512 3 

8.79 x 10 6 

0.5 


virialised matter contributes significantly to the non-linear 
clustering signal. 

Here, we examine the contributions of halo and non¬ 
halo mass to the matter power spectrum with the use of a 
set of N-body simulations. This paper is organised as follows. 
In l[2]we describe our simulations and the employed power 
spectrum estimator. In ([3] we investigate the contribution to 
the redshift zero matter power spectrum of haloes that are 
defined analogous to the typical halo model approach. We 
start by looking at the fraction of mass that is in haloes as 
a function of minimum halo mass and compare to analytic 
results in HU Next, in HI we examine the contributions 
of matter in regions with lower overdensities and outside of 
haloes, as a function of Fourier scale. We also examine what 
changes when we expand the haloes to include all matter 
associated to Friends-of-Friends (FoF) groups. I 11 43.2.11 we 
make predictions for the contribution of halo matter to the 
power spectrum as a function of both scale and minimum 
halo mass, which can serve as a test for halo models aimed 
at reproducing the clustering of dark matter. Finally, we 
summarise our findings in JU 


2 METHOD 
2.1 Simulations 


We base our anal ysis on a set of dark matter only runs 
from the OWLS (ISchave et all l20ld ~) and cosmo-OWLS 
iLe Brun et al.l [2014 1 projects. The simulations were run 
with a modified version of GADGET ill, the smoothed-particle 
h ydro dynamics (SPH) code last described in ISpringel 
(I2005T) . The cosmological parameters are derived from the 
Wilki nson Microwave Aniso tropy Probe (WMAP) 7-year re¬ 
sults ([Komatsu et al.|[201ll ~). and given by {9 m , 9b, 9 a, ers, 


n s , h} = {0.272, 0.0455, 0.728, 0.81, 0.967, 0.704}. 

We genera te_initial conditions assuming the 

lEisenstein fc Hul (11998)4 transfer function. Prior to im¬ 
posing the linear input spectrum, the parti cles are set u p 
in an initially glass-like state, as described in 1 Whit j (| 19941 ) . 
The partic l es are then evolved to redshift 2 = 127 using the 
IZel’dovichl (jl970l l approximation. 

The relevant parameters of the simulations we employ 
here are listed in Table [T| The simulation volumes range 
from 25/i -1 Mpc to 400/iT 1 Mpc. The mass resolution im¬ 
proves by a factor of 8 with each step, corresponding to 
an improvement of the spatial resolution by a factor of 2, 
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from the largest down to the smallest volume. The gravita¬ 
tional forces are softened on a comoving scale of 1/25 of the 
initial mean inter-particle spacing, L/N, but the softening 
length is limited to a resolution-dependent maximum phys¬ 
ical scale which is reached at 2 = 2.91. As we will demon¬ 
strate, by combining these simulations, we can accurately 
determine the matter power spectrum from linear scales up 
to k ~ 100 h Mpc -1 . 

2.2 Power spectrum calculation 

The matter power spectrum is a measure of the amount of 
structure that has formed on a given Fourier scale k, re¬ 
lated to a physical scale A through k = 2 -k/\. It is defined 
through the Fourier transform of the density contrast, Uk- 
We will present our results in terms of the dimensionless 
power spectrum, defined in the usual way: 



with V the volume of the simulation under consideration. 
As all particles have the same mass, the shot noise is simply 
equal to < |<5k| 2 >k,shot= 1/IVp, with N p the number of parti¬ 
cles in the simulation. All power spectra presented here have 
had shot noise subtracted to obtain more accurate results on 
small scales. 

We calculate the matter power sp ectrum using the pub - 
licly available f90 package powmes (IColombi et al.l 120091 ). 
The advantages of powmes stem from the use of the Fourier- 
Taylor transform, which allows analytical control of the bi¬ 
ases introduced, and the use of foldings of the particle dis¬ 
tribution, which allow the dynamic range to be extended 
to arbitrarily high wave numbers while keeping the statisti¬ 
cal errors b ounded. For a full des cription of these methods 
we re fer to IColombi et al.1 |2009). As in Ivan Daalen et all 
(120111) . we set the grid parameter to G = 256 and use a 
folding parameter F = 7 for the two smallest volumes. To 
calculate the power spectrum down to similar scales for the 
200 and 400 h -1 Mpc boxes, we set F equal to 8 and 9, 
respectively. Our results are insensitive to this choice of pa¬ 
rameters. 

Both box size and resolution effects lead to an under¬ 
estimation of the power - at least on scales where a suf¬ 
ficient number of modes is available so that the effects of 
mode discreteness can be ignored (fc > 8n/L) - while all 
simulations show excellent agreement on scales where they 
overlap (see Figure [21) ■ In order to cover the dynamic range 
from k = 0.01 h Mpc -1 to 100/iMpc -1 , we therefore com¬ 
bine the power spectra of different simulations by always 
taking the largest value of A 2 (fc) at each k. In the case of 
the full power spectrum, i.e. the power spectrum of all mat¬ 
ter, we take the combined power spectrum to be the one 
predicted by linear theory up to k = 0.12/iMpc -1 , where 
the power starts to become non-linear. While the largest 
boxes show excellent agreement with the linear power spec¬ 
trum on these scales, we wish to avoid box size effects as 
much as possible. For k > 0.12 /iMpc -1 - or, in the case of 
power spectra of subsets, for the smallest k- value available 
- we individually average each power spectrum over each of 
25 bins ki spaced equally in Fourier space and assign the 
combined power spectrum the largest A 2 {ki) derived in this 
manner between all simulations. 


We combine the power spectra of selections of particles 
(e.g. all particles that reside in haloes above a certain mass) 
in a similar way, but without including the linear theory 
power spectrum. 

Finally, we note that we take the contribution of halo 
matter to the power spectrum to be the auto-correlation of 
halo matter only (i.e. we do not examine the cross terms of 
halo and non-halo matter). 

2.3 Halo particle selection 

In the halo model approach, haloes are commonly defined 
through a spherical overdensity criterion, usually relative to 
the mean density of the Universe. In order to investigate the 
contribution of such haloes to the matter power spectrum, 
we define our haloes consistently. 

Overdense regions are identified in our simulations us¬ 
ing the Friends-of-Friends algorithm (with linking length 
0.2 times the mean interparticle distance), combined with a 
spherical o verdensity finde r, as implemented in the subfind 
algorithm (|Springel et alJl2001D . The centre of each region 
is taken to be the minimum of the gravitational potential. 
We define a halo as a spherical region with an internal mass 
overdensity of 200 x fl m p C rit, where p C rit is the critical density 
of the Universe. These haloes therefore have a mass equal 
to: 

47T q 

A/ 2 OO — Af200,mean = 200 X flmPcrit R 2 OO > (2) 

where R200 = i? 200 ,mean is the radius of the region. In the 
remainder of the paper, we will define halo particles as any 
particle with a distance R < R200 from any halo centre. 
All other particles are treated as non-halo particles, irre¬ 
spective of their possible FoF group membership, or having 
been identified as part of a bound subhalo by subfind. 

While we focus on halo matter as defined through R200 > 
we will also briefly discuss the contribution of halo matter to 
the power spectrum for other overdensity regions and halo 
definitions (i.e. R500, R2500 , fl^oo.crit and Friends-of-Friends) 
during the course of the paper. 

3 RESULTS 

3.1 Fractional mass in haloes 

We first examine the fraction of the mass that resides in 
haloes, fh- As in each simulation there is a lower limit to 
the masses of haloes that we can reliably resolve, we com¬ 
pute fh as a function of the minimum mass of the included 
haloes. Knowing the minimum resolved masses also allows 
us to estimate over which halo mass range we can probe the 
contribution of halo particles to the power spectrum in each 
simulation. 

The results for fh are shown in Figure [l] Different 
colours are used for each of our four different simulations, 
as indicated in the legend. Vertical dotted lines denote the 
masses corresponding to 100 particles. Below this limit the 
fraction of mass in haloes flattens off, indicating that such 
low-mass haloes are unresolved. A thick dashed line shows 
the result of combining the mass fractions of all four simula¬ 
tions for M min > 10 9 h -1 Mq, through /h.comb = max(/ h ,;), 
which we consider our best estimate for the true fh- The 
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| og,o( M m;„/[ M 0 / h ]) 

Figure 1. Top panel: The cumulative fraction of mass inside 
haloes, /h, as a function of minimum halo mass, for different col¬ 
lisionless simulations as indicated in the legend. The resolution 
limit, defined as the mass of haloes containing 100 particles, is 
shown as a vertical dotted line for each simulation. Below this 
limit, the fraction of mass in haloes is underestimated. For the 
two highest-resolution simulations ( L050 and L025) these frac¬ 
tions are also significantly underestimated at high masses, as such 
haloes are under-represented in these small volumes. Between the 
limits imposed by resolution and box size effects, the simulations 
are in excellent agreement, and show that the fraction of mass 
in haloes is ~ 52% for M 200 > 10 9 h~ 1 Mq. The black dashed 
line shows the combined result, taking the maximum fraction of 
mass in haloes between the diffe rent simulations at every mass. 
We also show predictions for the Ti nker et alj J2008T ) mass func¬ 
tion as a black dot-dashed line (see main text). Bottom panel: 
The fraction of this combined function, fh/fh combi predicted by 
each simulation. 


bottom panel shows the ratio of /h of each simulation to 
this combined fraction. 

At the massive end, the high-resolution but low-volume 
L025 and L050 simulations significantly underestimate /h- 
This is most clearly seen in the bottom panel: for L025 
the mass fraction in haloes is significantly underestimated 
for halo masses M 200 > 10 11 (i _1 Mq, while for the L050 
box this happens for M 200 > 10 12 /i -1 Mq. These values 
correspond to the masses above which the halo mass func¬ 
tions are underestimated for these simulations (not shown). 
The fluctuations seen in the bottom panel for L200 where 
A /200 > 10 14 h -1 Mq are due to the rarity of such mas¬ 
sive haloes, but as the fraction of the mass residing in such 
haloes is < 10% this does not impact our conclusions. All 
simulations in which haloes at a certain A/ m i n are both well- 
resolved and well-represented show excellent agreement for 
fh{M > M m in). 

The fraction of mass in haloes increases with decreas¬ 
ing halo mass. Only ~ 19% of matter is found in groups and 
clusters (A/ m i n > 10 13 h _1 Mq), which increases to ~ 30% 
for haloes with masses greater than that of the Milky Way 
(A/ m in > 10 12 hT 1 Mq). But even at the lowest resolved 
mass of roughly 10 9 h~ 1 Mq, the fraction of mass in haloes is 


still barely more than 50%. We therefore expect a significant 
contribution from particles in haloes with M < 10 9 / 2, -1 Mq 
- and possibly from dark matter particles that do not reside 
in haloes of any mass - to the matter power spectrum on 
large scales, which we calculate in the next section. 

For comparison, the top panel of Figure [T] also shows 
predictions f or the fraction of mass in haloes from the 
iTinker et al.l ( 2008 ) M 200 halo mass function. Using the nor¬ 
malized halo mass function fit provided by these authors, 
we have calculated fh(M > A/ m i n ) under the standard halo 
model assumption that all mass resides in haloes (i.e. the 
fit converges to unity when all halo masses are included). 
The results are shown by the black dot-dashed line. Up to 
Mm in ~ 10 12 /i _1 Mq the relative difference between the 
ITinker et al.l (120081) prediction and our combined result is 
constant at about 10% before decreasing at higher masses. 
One possible reason for this discrepancy is that we count 
matter in regions where haloes overlap only once, whereas 
double counting is possible when integrating the mass func¬ 
tion. However, we have checked that the mass residing in 
overlap regions in our simulation is always < 1.7%, with 
the largest overlap fraction being found for the most mas¬ 
sive haloes. The < 10% differences found for /h are there¬ 
fore likely due to the non-universality o f the halo mass 
function at this level of precision (e.g. ITinker et al .1 EfOOsl . 
iMurrav. Power fc Robothamll20i3 ). or perhaps cosmic vari¬ 
ance. 

Whether or not fh(M < A/ m i n ) converges to unity in re¬ 
ality depends on the nature of dark matter, but for perfectly 
cold dark matter the expectation is that it should. We note, 
though, that the convergence of /h with mass is extremely 
slow. Taking the [Tinker et alj ( 2008 ) fit as an example, the 
fraction of mass in haloes with A/200 > 10 9 h” 1 Mq is 0.56, 
and this number is still only 0.76 for A/200 > 1 h _1 Mq. Even 
for A/200 > 1kg, fh(M < Mmin) ~ 0.88. The prediction that 
a significant amount of dark matter is in ultra-small haloes 
means that the line between halo matter and truly smoothly 
distributed matter is vague - but as we will show, it is in at 
least some cases unnecessary to make the distinction. 

Since fh(M > M m in ) continues to rise down to mass 
scales that are unresolved by our simulations, we expect to 
underestimate the total contribution of matter in haloes to 
the power spectrum. However, as we will see in f )3.2.1l this 
depends on the spatial scale considered. There exists a range 
in Fourier space where the fraction of power from halo par¬ 
ticles converges to values below unity, and the contribution 
from haloes with masses A/200 10 11 h ~ 1 Mq is negligible. 

On scales where this does not hold we can still constrain the 
contribution from haloes above a certain mass. 

In the remainder of the paper, we will only consider 
particles residing in haloes with A/200 > 10 9 h _1 Mq to be 
halo particles, as this corresponds roughly to the smallest 
haloes we can resolve. 

3.2 Halo contribution to the power spectrum 

We first show the full dimensionless matter power spec¬ 
trum, i.e. using all particles, in Figure [2] Here each simula¬ 
tion is shown using a different colour, and it is immediately 
clear that no single one is converged over the full range of 
wavenumbers. The linear theor y power spectrum, as gener- 
ated by the f90 package CAMB (|Lewis. Challinor fc Lasenbvt 
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A [Mpc/h] 



0.1 1.0 10.0 100.0 1000.0 
k [h/Mpc] 

Figure 2. Top panel: The dimensionless power spectrum de¬ 
rived from each simulation, along with the linear power spec¬ 
trum (long-dashed purple line) and the combined power spec¬ 
trum (dashed black line). While the L025 and L050 simulations 
significantly underestimate the power on large scales due to miss¬ 
ing modes, their high resolution allows us to accurately extend 
the power spectrum obtained using the larger volumes up to 
k ~ 100/iMpc -1 . The erratic behaviour seen for low-resolution 
simulations at large k is due to shot noise subtraction. Bottom 
panel: The fraction of power relative to the combined power spec¬ 
trum for each simulation (as well as for the linear theory predic¬ 
tion). For k < 20hMpc —1 , multiple simulations show the same 
result, indicating convergence on these scales. 


l200d. version January 2010), is shown as the long dashed 
purple line. Simulations L400 and L200 show good agree¬ 
ment with the linear power spectrum on scales where non¬ 
linear evolution is negligible (fc < 0.12 hMpc -1 ) and a 
sufficient number of modes is available (k > 0.04 and 
0.08 hMpc -1 respectively, roughly corresponding to A = 
0.4 L), while L050 and L025 show severe box size effects 
due to their lack of large-scale modes. These box size ef¬ 
fects become negligible only for k > 10 and k > 40/iMpc -1 , 
respectively. 

Due to their finite resolution, all simulations underes¬ 
timate the power on small scales. Note that shot noise was 
subtracted from all power spectra shown here, which ex¬ 
plains the erratic behaviour of the power spectra on the 
smallest scales. The underestimation of small-scale power 
becomes significant already on scales corresponding to ~ 100 
softening lengths. However, for every wave number k < 
100/iMpc -1 , there is at least one simulation for which nei¬ 
ther box size nor resolution leads to an underestimation of 
the power at the > 1% level. We therefore combine the 
different power spectra as described in m\ to obtain the 
combined power spectrum, P com b = max(P;), shown as the 
dashed black line. 

The bottom panel of Figure [2] shows the fraction of 
power predicted by each simulation, as well as the fraction 
predicted by linear theory, relative to the combined power 


A [Mpc/h] 

100.0 10.0 1.0 0.1 



Figure 3. The combined power spectrum for different sets of 
particles: R < R 200 (halo particles), R < R500, R < R 2500 and 
R > -R200 (non-halo particles). Only haloes with A/200 > M m ; n = 
10® hr 1 Mg were considered in the cuts made, which in total con¬ 
tain about 52% of all dark matter. The halo particles easily dom¬ 
inate the power on small scales. On linear scales, the contribution 
of halo and non-halo particles is roughly equal, which is expected 
as the relative power on these scales depends mostly on the mass 
in each component. The cross-terms between the halo and non¬ 
halo particles (not shown) account for about half the power in 
the linear and mildly non-linear regimes (fc < 0.2 ftMpc -1 ). Note 
that the horizontal range has been shortened relative to Figure [2] 


spectrum. By construction, this fraction is bounded to unity 
on non-linear scales. Note that on scales k < 20 h Mpc^ 1 , the 
fractions of multiple simulations are within a few percent 
of unity, indicating convergence on these scales. For smaller 
scales, however, convergence is uncertain, although based on 
the results for larger scales we expect our combined power 
spectrum to be accurate to ~ 1% up to k ~ 100 hMpc^ 1 . 

Next, we repeat this procedure for halo and non¬ 
halo particles. We also consider particles within the -R 500 
and -R 2500 overdensity regions, defined analogously to R200, 
which probe the inner parts of haloes. As we cannot reli¬ 
ably resolve haloes with less than about 100 particles in any 
simulation, we only consider the contribution of haloes with 
masses M200 > M m i n = max [10 9 h _1 Mg, 100 mam] here, 
referring to matter in lower-mass haloes as non-halo parti¬ 
cles. The results are shown in Figure[3] Note that for clarity 
only the combined power spectra are shown, and that the 
horizontal range has been shortened with respect to Fig¬ 
ure [5] only showing the range of scales for which we can re¬ 
liably determine the power spectrum. Furthermore, we note 
that the power spectrum of each component (e.g. particles 
within -R 500 , particles outside of R200) was calculated inde¬ 
pendent of the total matter and not renormalized, in order 
to facilitate a direct comparison of the power contained in 
each as a function of scale. 

The contribution from halo particles strongly dominates 
the power on small scales. The halo contribution is in turn 
dominated by the very inner regions of haloes, at least on 
scales smaller than the size of these regions. However, to- 
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X [Mpc/h] 

100.0 10.0 1.0 0.1 



Figure 4. The fraction of power within haloes with masses 
M 200 > 10 9 h~ 1 Mq as a function of wavenumber. The dashed 
black line shows the combined power spectrum derived from the 
smoothed power spectra of the four simulations employed in this 
paper, each of which is shown as well. The halo contribution 
rises rapidly down to A ~ 2/i _1 Mpc, peaking at ~ 95% for 
k « 20/iMpc —1 (A ~ 300 h~ 1 kpc) and remaining roughly con¬ 
stant for larger k. On smaller scales the power spectrum is dom¬ 
inated by increasingly smaller haloes, while on the largest scales 
the contribution of haloes to the power spectrum depends mainly 
on the total mass fraction in haloes. 


X [Mpc/h] 

100.0 10.0 1.0 0.1 



Figure 5. As Figure |4] but now for all mass inside FoF groups 
with M 200 > 10 9 h~ 1 Mq (we select on M 200 in order to keep our 
halo sample identical to the one used in Fig. 0 . While the scale de¬ 
pendence is very similar (i.e. a rapid rise down to A ~ 1 h~ 1 Mpc 
and roughly constant on smaller scales), the contribution to the 
power spectrum is higher than for the R 200 overdensity regions 
(shown as a solid grey line) on any scale. The contribution of halo 
matter to the power spectrum is increased by 5 — 10% on most 
scales relative to the results of FigiireQ] and matter in FoF groups 
accounts for essentially all power on scales k > 3ftMpc _1 . This 
implies that the R 200 overdensity regions do not fully capture the 
halo. The grey dashed line shows the combined result if R 200 ,crit 
is used. 


wards larger scales this contribution diminishes, and for 
k < 0.4 h Mpc -1 less than half of the total power is pro¬ 
vided by matter in haloes alone. On large scales the sig¬ 
nificant fraction of the mass that does not reside in haloes 
with M 200 > 10 9 /2, -1 Mq becomes more important, its con¬ 
tribution to the power spectrum increasing to about 20%, 
almost equalling the contribution of halo matter on linear 
scales. The remaining ~ 40% of the total matter power on 
large scales is therefore contributed by the cross-terms of 
halo and non-halo matter (not shown here). 

Note that on the scales shown here, only L400 and L200 
contribute to the combined power spectrum of non-halo par¬ 
ticles. Nonetheless, as Figure [2] shows that these two simu¬ 
lations are in excellent agreement for 0.4 < k < 10 h Mpc -1 
even though the mass resolution is eight times worse for 
L400, we have no reason to believe that this component 
would change significantly on non-linear scales if lower-mass 
haloes were resolved. On linear scales, however, the contri¬ 
bution of halo matter is mostly determined by the fraction 
of mass in haloes, which does depend on the minimum halo 
mass resolved. We will return to this point in 43.2.11 

We investigate the contribution of halo matter in more 
detail in Figure [4] which shows the ratio of the power spec¬ 
trum of matter within R 200 of haloes with masses M 200 > 
10 9 h -1 Mq to the power spectrum of all matter. The black 
dashed line shows the ratio of the combined power spec¬ 
tra, obtained from the smoothed power spectra of all four 
simulations shown here as described in 42.21 relative to the 
combined total power spectrum (black line in Figure [3]). The 


solid lines show the relative contributions of halo matter sep¬ 
arately for each simulation. 

The contribution of halo matter to the total power 
increases with decreasing physical scale. On large (linear) 
scales, the contribution from haloes seems to converge to ~ 
30%, in good agreement with /h (M > 10 9 /i _1 Mq) 2 « 0.27. 
This is expected, as the contribution of any subset of mat¬ 
ter to the power spectrum on sufficiently large scales should 
scale only with (the square of) the fraction of mass con¬ 
tained in such a subset. However, as the fraction of power 
in haloes on large scales is fully determined by L200 and 
L400, with both predicting roughly the same fraction as can 
be seen in Figure [4] while the fraction of mass in haloes 
M > 10 9 hT 1 Mq is only accurately measured for L025, this 
correspondence is actually surprising. 

On non-linear scales the ratio of the power from halo 
matter to the total power increases rapidly down to phys¬ 
ical scales of A ~ 2ft -1 Mpc (k ~ 3/iMpc -1 ), reaching at 
most 95%, before slowly levelling off towards smaller scales. 
Note that the combined results are fully determined by L050 
around k « 20 h Mpc -1 , where we are unable to show con¬ 
vergence because the resolution of L200 is too low and the 
volume of L025 is too small. However, we will show in 43.2.11 
that on these scales little would change if higher-resolution 
simulations were available. 

While L400 and L200 are in good agreement for 0.2 < 
k < 10 h Mpc -1 , on sub-Mpc scales the contribution of halo 
matter to the total matter power spectrum starts to show 
a strong dependence on resolution. On these scales fluctua- 
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tions within the same halo (i.e. the 1-halo term in halo model 
terminology) dominate the power spectrum, so naturally the 
contribution to the power will be underestimated on scales 
A < -R 200 ,min, where f? 200 ,min is the virial radius of a halo 
with the minimum resolved mass, M m i n , in that particular 
simulation. In practice, the power is already significantly 
underestimated on larger scales. Fortunately, the combina¬ 
tion of simulations chosen here still allows us to probe the 
contribution of halo matter up to fc max ~ 100 ftMpc -1 . 

As the power on sub-Mpc scales is predicted to be dom¬ 
inated by the 1-lialo term (which we also observe in our sim¬ 
ulations on resolved mass scales), adding lower-mass haloes 
than those resolved here is expected to have a negligible im¬ 
pact on the measured contribution of halo matter on scales 
fcmax < 27r/I?200,min. Therefore, 5 — 7% of small-scale power 
is unaccounted for by halo particles, regardless of resolu¬ 
tion effects. Note that it is possible that the 2-halo corre¬ 
lations between unresolved and other haloes are responsi¬ 
ble for making up the deficit. However, to explain our re¬ 
sults these unresolved haloes would have to cluster directly 
around the resolved haloes, from which a picture arises that 
is essentially the same as viewing “smooth” halo matter as 
being made up entirely of tiny haloes. We find that it is the 
cross-term between halo matter and matter just outside the 
R 200 regions that makes up the deficit. 

To demonstrate that this is indeed the case, we calcu¬ 
late the contribution of matter in FoF groups to the total 
power spectrum, with a mass limit of M 200 > 10 9 ft -1 Mq. 
Note that we select on M 200 in order to keep our halo sample 
identical. The results are shown in Figure [5] The combined 
result of Figure [4] is shown as a solid grey line to aid the 
comparison. While the scale-dependence of the halo contri¬ 
bution for FoF groups is similar to that shown for the R 200 
regions, the halo contribution is significantly larger on all 
scales (with the exception of scales k ~ 2 ftMpc -1 ), and is 
essentially 100% for k > 3 ftMpc -1 . This Fourier scale cor¬ 
responds to the virial radius of the largest clusters in the 
simulation. 

On scales k < 0.3 ftMpc -1 , the contribution of matter 
in FoF groups to the clustering signal is consistently ~ 20% 
higher than that of matter in R 200 haloes. Interestingly, the 
fraction of mass in FoF groups is only about 4% higher than 
that in R 200 haloes (not shown). This implies that the ob¬ 
served increase in the contribution of halo matter to the 
power spectrum when using the FoF instead of the R -200 re¬ 
gion is not only due to the addition of mass, but mainly due 
to the addition of clustered material. 

Finally, we also show the results if I?2oo,crit is used in¬ 
stead (still with M 200 > lO 9 ft -1 M 0 ), as a dashed grey 
line in Figure [5] As such an overdensity criterion picks out 
smaller regions than R 200 , containing less mass, the contri¬ 
bution of halo matter to the power spectrum is also smaller, 
especially on large scales. On sub-Mpc scales, however, the 
differences are small, with the contribution to the power 
spectrum of halo matter peaking at 94%. 

We conclude that what region is chosen to represent a 
halo has a large impact on the contribution of haloes to the 
matter power spectrum, in a scale-dependent way. In what 
follows, we will continue to define haloes using the mean 
overdensity criterion, as this is typically used in the halo 
model approach. 


3.2.1 Mass dependence 

To see which halo masses contribute most to the matter 
power spectrum as a function of scale, while simultaneously 
examining the dependence of our results on the mass of the 
lowest resolved halo, we turn to Figure [U] Each panel corre¬ 
sponds to a different simulation and each curve to a different 
minimum halo mass. The halo contributions are shown rel¬ 
ative to the combined power spectrum of all matter (black 
line in Figure [3]). 

The legend shows the minimum halo mass, 
fog lo (M min /[M 0 /ft,]), that corresponds to each curve. 
Note that the minimum masses differ for each simulation, 
because they start at 100 particles. For each simulation the 
minimum halo masses of the different curves are half a dex 
apart. 

Grey regions indicate the approximate scales on which 
the full matter power spectrum of the simulation is not con¬ 
verged to ~ 1% with respect to the combined one. While 
this gives an indication of which scales to trust, note that 
the relative contribution of each halo mass can be converged 
on a different range of scales. 

Finally, the bottom half of each panel shows the differ¬ 
ence between consecutive curves, i.e. the relative contribu¬ 
tion added by decreasing the minimum halo mass by half a 
dex. Here /a,* = A| 0 o,i/A a n. As we will show shortly, while 
the relative contributions of haloes of a certain mass shown 
in the bottom halves of the panels can be compared between 
different simulations, the same does not hold for the abso¬ 
lute contributions, as box size effects play an important role 
on a large range of scales. 

As Figure [6] shows, all haloes provide a significant 
contribution to the power on the largest scales (k < 
0.2 ftMpc -1 ). It is clear that we cannot claim convergence 
on these scales. On sub-Mpc scales, low-mass haloes become 
increasingly important as one moves to larger values of k. 
This means that, as expected, one needs to resolve smaller 
haloes to obtain convergence on smaller scales. Interestingly, 
there is a region in between these two regimes where we do 
see convergence with decreasing halo mass. This is most eas¬ 
ily seen in the bottom half of each panel, where the conver¬ 
gence manifests itself as a sharp drop to very low values for 
low halo masses, indicating that the addition of lower-mass 
haloes has a negligible effect on the power spectrum. 

Focusing on the L200 simulation (top right panel), 
we see that it is converged with minimum halo mass on 
scales k ~ 4 ft. Mpc -1 , roughly where the contribution from 
halo matter plateaus. On these scales, haloes with M 200 < 
10 11 ft -1 M 0 provide a negligible contribution to the power 
spectrum. This is confirmed by looking at the bottom half 
of the bottom left panel, showing the results for the L050 
simulation: haloes below ~ 10 11 ft -1 M 0 do not measur¬ 
ably impact the power for k ~ 4 ft Mpc -1 . As we can also 
see from this panel, the range of scales on which conver¬ 
gence is obtained widens as the minimum resolved halo mass 
decreases. This convergence, coupled with the result that 
A|oo(^)/ Atot (fc) < 1, shows that on some scales not all 
power comes from matter inside f? 200 - 

One might worry that the effects of box size skew these 
results. To test this, we turn to Figure [7] which is identi¬ 
cal to the top right panel of Figure [6] but with the results 
of LI00 superimposed as dashed lines. The LI00 simula- 
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Figure 6. Comparison of the contribution of haloes above a given mass to the matter power spectrum relative to the total combined 
power spectrum of all simulations, for L400 (top left), L200 (top right), L050 (bottom left) and L025 (bottom right). The legend shows 
the minimum halo mass log 10 (M m i n /[MQ//i]). Note that lines of the same colour do not correspond to the same minimum halo mass 
in the four panels, as the binning is based on the minimum resolved halo mass (see text). The grey regions indicate where box size or 
resolution effects are ->1% for the full power spectrum; the relative contribution of specific halo masses may be converged on a different 
range, as can be seen by comparing the panels. The bottom half of each panel shows the difference between consecutive curves in the 
top panel, i.e. the relative contribution added by decreasing the minimum halo mass by half a dex. Several panels show convergence of 
the fractional power with M min to values smaller than unity for wavenumbers that are not greyed out, indicating that the conclusion 
that matter outside haloes contributes significantly to the power spectrum is robust. 
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k [h/Mpc] 

Figure 7. As the top-right panel of Figure [6] but with the re¬ 
sults for LI00 added as dashed lines for the same minimum halo 
masses, showing the effects of box size at fixed resolution. Due 
to the missing large-scale modes in LI00, the large-scale con¬ 
tribution is underestimated. Additionally, high-mass haloes are 
under-represented and the role of low-mass haloes on small scales 
is overestimated. However, the relative contributions of haloes of 
a certain mass shown in the bottom half of the figure are in ex¬ 
cellent agreement for all but the highest mass bin. 


tion has 512 3 particles and therefore the same resolution 
as L200 , but in an 8 x smaller volume. Comparing the two 
simulations therefore shows the effects of box size at fixed 
resolution. On large scales and for high-mass haloes, the 
contribution of halo matter is underestimated in LI00, rel¬ 
ative to L200. This is expected, as large-scale modes are 
missed in the smaller boxes and massive haloes are under¬ 
represented. Meanwhile, the contribution of low-mass haloes 
on small scales tends to be overestimated, even though the 
resolution is identical. Interestingly, there are mass and spa¬ 
tial scales where the simulations are in near perfect agree¬ 
ment, such as for a minimum halo mass of 10 13 ' 75 h^ 1 Mq 
and k > 3hMpc -1 . Most important, however, is that the 
contributions from haloes above a given halo mass shown in 
the bottom half of the panel are in excellent agreement for 
the two simulations over the entire range of scales, except 
for the highest mass bins (which are under-represented in 
L100 ) and the largest modes. This shows that we can still 
derive the correct contribution of haloes within a certain 
mass range, and investigate whether we are converged with 
mass on a certain scale, even when box size effects play a 
role. For example, while the simulations shown in Figure [6] 
do not predict exactly the same contribution to the power 


spectrum for haloes with M 200 > 10 11 ft -1 Mq, due to the 
limited box size of L050 and L025, it is clear that the results 
for the relative importance of such haloes are converged. In 
Appendix 0 we investigate the effects of box size and reso¬ 
lution further, and conclude that the contribution of haloes 
in a certain mass range to the power spectrum is indeed 
converged with box size over an interestingly large range of 
scales, while resolution only affects the conclusions on small 
scales. 

Figure [6] makes other interesting predictions as well. 
For example, comparing again the results for L200 and 
L050, we see that both simulations agree that most of the 
power at k ~ 10/iMpc^ 1 comes from haloes with masses 
M 200 > 10 13 / 2 T 1 Mq. These group and cluster-scale haloes 
remain the dominant contributors on somewhat larger scales 
as well, their contribution peaking around k = 2 — 3 hMpc -1 
before gradually falling off. Note that haloes with M 200 > 
10 13 h _1 Mq only account for about 19% of the total mass 
(see Figure |T}). This provides an interesting test for halo 
models, where the contribution of haloes above a certain 
mass is strongly dependent on the adopted concentration- 
mass relation. 


4 SUMMARY & CONCLUSIONS 

In this work we investigated the contribution of haloes to 
the matter power spectrum as a function of both scale and 
halo mass, thus testing the assumption typically made by 
halo-based models that all matter resides in spherical haloes. 
To do so, we combined a set of cosmological N-body simu¬ 
lations to calculate the contributions of different spherical 
overdensity regions, FoF groups and matter outside haloes 
to the power spectrum, paying careful attention to the con¬ 
vergence with both numerical resolution and the size of the 
simulation volume. As convergence with mass is generally 
very slow, any claims about the role of haloes or the mass 
contained in them need to be quoted together with a min¬ 
imum halo mass in order to have a meaningful interpreta¬ 
tion. We note that when we refer to the contribution of halo 
matter to the matter power spectrum, we consider only the 
auto-correlation of halo matter, ignoring the cross terms of 
halo and non-halo matter. 

Our findings can be summarised as follows: 

• On scales k < 1/iMpc -1 , haloes - defined as spheri¬ 
cal regions with an enclosed overdensity of 200 times the 
mean matter density in the Universe - with masses M 200 < 
10 9 ' 5 ft -1 Mq, which are not resolved here, may contribute 
significantly to the matter power spectrum. For 2 < k < 
60/iMpc^ 1 , our simulations suggest their contribution to 
be < 1%. A range of scales around k ~ 4/iMpc -1 exists for 
which the contribution of halo matter to the power spectrum 
appears to be fully converged with decreasing halo mass. 

• Matter within -R 200 alone cannot account for all power 
f° r 1 tS k 100/rMpc^ 1 . Its relative contribution increases 
with increasing Fourier scale, peaking at ~ 95% around 
k = 20/iMpc _1 . On smaller scales, its contribution remains 
roughly constant. 

• When Jfeoo.crit is used to define haloes instead of the 
fiducial R 200 , the contribution of haloes to the power spec¬ 
trum decreases significantly on all scales. 
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• Matter just outside the -R 200 , but identified as part 
of FoF groups, provides an important contribution to the 
power spectrum. Taken together, matter in FoF groups with 
Mfof > 10 9 h Mg accounts for essentially all power for 
3 < k < 100 ft Mpc -1 . Switching from R 200 to FoF haloes 
increases the contribution of halo matter on any scale probed 
here by 5 — 15%. 

• For 2 < k < 10 hMpc -1 , haloes below ~ 10 11 fi _1 Mg 
provide a negligible contribution to the power spectrum. The 
dominant contribution on these scales is provided by haloes 
with masses M 200 > 10 13 ' 5 ft -1 Mq , even though such haloes 
account for only ~ 13% of the total mass. 

As we have demonstrated, the halo model assumption 
that all matter resides in (spherical overdensity) haloes may 
have significant consequences for the predictions of the mat¬ 
ter power spectrum. Specifically, we expect such an approach 
to overestimate the contribution of haloes to the power on 
small scales (k > 1 ftMpc -1 ), mainly because it ignores the 
contribution of matter just outside R 200 to the power spec- 
trum[f| While defining haloes to be larger regions similar 
to FoF groups mitigates the small-scale power deficits, the 
fact that such regions are often non-virialised and typically 
non-spherical may lead to other problems. 

Clearly, the validity of the postulate that the clustering 
of matter is fully determined by matter in haloes is strongly 
dependent on the definition of a halo used - but it is hard 
to say what the “best” definition to use in this context is. 
For example, while haloes defined through i? 200 ,crit will be 
more compact and therefore have a smaller overlap fraction 
than R 200 or FoF haloes, their contribution to the power 
spectrum will be smaller for the same minimum halo mass. 
And while FoF groups contain nearly all the mass that is im¬ 
portant for clustering on small scales, the fact that they are 
not completely virialised, are non-spherical and have bound¬ 
aries tha^_do_jio^_correspond to a fixed mean overdensity 
(e.g. iMore et al.lboill ') prohibits their use in traditional halo 
based models. 
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APPENDIX A: ADDITIONAL CONVERGENCE 
TESTS 

Box size and resolution effects can have a large impact on the 
power spectra measured from simulations. With decreasing 
box size the number of large-scale modes decreases, lead¬ 
ing to an underestimation of the power. Massive haloes also 
become under-represented, which may lead to an underes¬ 
timation of the contribution of haloes above a certain mass 
to the power spectrum. Conversely, with increasing particle 
mass the minimum mass at which haloes can be reliably re¬ 
solved increases as well. Here, we investigate these effects 
independently by changing the box size at fixed resolution 
and vice versa. 

In Figure [All we show the effects of decreasing the box 
size while keeping the resolution fixed. We do this for two 


different resolutions: high resolution, based on L200, on the 
left-hand side, and low resolution, based on LJ^OO, on the 
right-hand side. For clarity the curves were smoothed by im¬ 
posing a minimum bin size in k of 0.05 dex and averaging the 
power over all modes that fall in each bin. The results for the 
largest boxes are shown with solid lines, the intermediate- 
size boxes with dashed lines, and the smallest boxes with 
dot-dashed lines. The top panels show the relative contribu¬ 
tions of haloes above a certain mass to the power spectrum, 
while the bottom panels show the relative contributions of 
haloes in mass ranges 0.5 dex wide. 

Looking at the top panels, we see that changing the box 
size can severely affect the derived contribution of haloes 
above a given mass. Not only does decreasing the box size 
lead to a large underestimation of both the power on large 
scales and of the contribution of massive haloes (especially 
for L050N256), but the power on small scales is simulta¬ 
neously overestimated in the smaller boxes. However, if we 
instead turn to the bottom panel we see that the contri¬ 
butions of haloes in a certain mass range are much better 
converged and over almost the entire range, with the ex¬ 
ception of the principal modes and the highest mass bin in 
both figures. The relative contribution of haloes in a certain 
mass range can therefore be used to investigate convergence, 
and is the preferred quantity to compare against halo model 
predictions. 

Finally, in Figure l~\~2l we decrease the resolution while 
keeping the box size fixed. We show results for simulations 
with a 200 h -1 Mpc box on the left-hand side, and for simu¬ 
lations with a 100 h -1 Mpc box on the right-hand side. Both 
panels show that decreasing the resolution only has a sig¬ 
nificant effect below some length scale and below some halo 
mass. For any halo mass limit a length scale exists down to 
which the power spectra of simulations with different reso¬ 
lutions show excellent agreement. 

This paper has been typeset from a TgX/ DT^X file prepared 
by the author. 
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Figure Al. Convergence test where the box size is changed at fixed resolution, for three high-resolution (left) and three low-resolution 
simulations (right). The different simulations are shown as solid (large box), dashed (intermediate box) and dot-dashed (small box) lines. 
The top panels show the relative contribution of haloes to the power spectrum for different minimum halo masses, shown in the legend 
as log-j^(M m in/[M 0 /h]). The grey regions indicate where box size or resolution effects are > 1% for the full power spectrum. The curves 
were smoothed to improve their visibility (see text). The bottom half of each panel shows the difference between consecutive curves in 
the top panel, i.e. the relative contribution added by decreasing the minimum halo mass by half a dex. The effects of decreasing the box 
size are especially apparent for L050N256 and for high halo masses, which are under-represented in the smaller boxes. As we saw in 
Figure |3 lowering the resolution while keeping the box size fixed decreases the relative contribution to the power spectrum for haloes 
above some mass at large scales, while increasing it at small scales. However, the relative contributions of each range in halo mass (shown 
in the bottom panels) agree very well in each simulation (with the exception of the principal modes and highest mass bin), even when 
the box size effects are large. 
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Figure A2. As Figure lAll but now the resolution is changed at fixed box size, for simulations with box sizes of 200 A 1 Mpc (left) and 
100 X -1 Mpc (right). Note that we only show lines at masses where haloes are resolved with at least 40 particles. Lowering the resolution 
causes the power on small scales and in low-mass haloes to be underestimated, which affects both the relative contribution of haloes 
above a certain mass (top panels) and in a certain range of halo masses (bottom panels). However, for any halo mass limit a length scale 
exists down to which the power spectra of simulations with different resolutions show excellent agreement. 
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